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Abstract 

Rate processes with dynamical disorder are investigated within a simple framework provided 
by unidirectional electron transfer (ET) with fluctuating transfer rate. The rate fluctuations are 
assumed to be described by a non-Markovian stochastic jump process which reflects conformational 
dynamics of an electron transferring donor-acceptor molecular complex. A tractable analytical 
expression is obtained for the relaxation of the donor population (in the Laplace-transformed 
time domain) averaged over the stationary conformational fluctuations. The corresponding mean 
transfer time is also obtained in an analytical form. The case of two-state fluctuations is studied 
in detail for a model incorporating substate diffusion within one of the conformations. It is shown 
that an increase of the conformational diffusion time results in a gradual transition from the regime 
of fast modulation characterized by the averaged ET rate to the regime of quasi-static disorder. 
This transition occurs at the conformational mean residence time intervals fixed much less than 
the inverse of the corresponding ET rates. An explanation of this paradoxical effect is provided. 
Moreover, its presence is also manifested for the simplest, exactly solvable non-Markovian model 
with a biexponential distribution of the residence times in one of the conformations. The nontrivial 
conditions for this phenomenon to occur are found. 

PACS numbers: 05.40.-a,82.20.-w,82.20.Uv 
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I. INTRODUCTION 



Rate processes with dynamical disorder are of prominent importance in many areas of 
physics and chemistry 1 . Consider for example a photo-induced electron transfer (ET) in a 
donor-acceptor ET complex like bacterial photosynthetic center, or various artificial donor- 
bridge-acceptor molecular systems^. Light excites the electron donating molecule from its 
ground electronic state, D, to the excited state, D*. From D* the electron is transferred 
to the acceptor molecule A, D*A — > D + A~, from which it returns back to the donor in its 
ground state. Such ET process can be repeated many times. If the corresponding energy 
differences are large, the backward reactions are negligible. In the Golden Rule approxi- 
mation, the nonadiabatic ET rate is determined primarily by a strong electron-vibrational 
coupling and a weak tunneling coupling between the diabatic electronic states. It is given 
by k[£\ = ^|V^[£]| 2 FC[£], where V[£] is the effective tunneling matrix element and FC[£] is 
the corresponding Franck-Condon factor^- Both V[£] and FC[£] can depend in addition on 
a conformational coordinate which reflects stochastic conformational dynamics of the 
ET complex. Such dynamics is slow on the time scale of vibrational relaxation and the rate 
follows instantaneously to the conformational changes. Under these assumptions, the 
ET kinetics is described by a simple stochastic relaxation equation for the decay of donor 
population 

p D {t) = -k[£{t)]PD{t) (1) 

with a fluctuating rate &[£(£)]. This presents a simplest example of rate processes with 
dynamical disorder 1 . 

Similar descriptions have been introduced in ET theory by Agmon and Hopfield^ and 
others 5 . In these studies, the conformational fluctuations were assumed to be continuous 
obeying a Smoluchowski equation description. In such a case, the Kramers diffusion theory 
approach^ provides a standard theoretical framework. The conformational fluctuations can, 
however, consist of discrete, abrupt changes. Then, a discrete state Markovian modeling 
seems more appropriat e) 8 ! 9 ! 10 ! 11 . It should be however remembered that the corresponding 
discrete conformational states present in reality some manifolds of many substates with a 
possibly complex internal structur e) 12 ! 1 ? and the relevant intrinsic dynamics can produce 
non-exponential, broad residence time distributions (RTDs) of the time intervals Tj spent by 
the molecule in the corresponding discrete macrostates j. It can be, for example, a power 
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law^. For this reason, an approximate discrete state description should generally be non- 
Markovian. Similar non-Markovian features, which are revealed by non-Poissonian statistics 
in the dynamics of single molecules^, came recently in the focus of attentio n 15 ^ 6 . 

The main goal of this work is to clarify the role of non-Markovian memory effects for 
the rate processes with dynamical disorder. This challenging problem is exemplified by a 
two-state model incorporating a conformational diffusion inside one of the macro conforma- 
tions. It is shown that the non-Markovian effects can be of profound importance leading 
to a paradoxical effect. Namely, the averaged ET kinetics can obey a quasi-static disorder 
description notwithstanding that the ET rate fluctuates rapidly on the time scale of aver- 
aged ET process. A resolution of this apparent paradox is provided. Its presence is also 
demonstrated in the simplest, exactly solvable non-Markovian model with a biexponential 
distribution of the residence time in one of the conformations. It can be thus of general 
importance. 

II. MODEL AND THEORY 

Let us consider the relaxation dynamics in Eq. ((TJ with £(t) being a discrete state non- 
Markovian jump process termed also the renewal process or continuous time random 
walk (CTRW) 18 i 19 . It is characterized by the RTDs t/> 3 -(t) and the matrix P^ of the transition 
probabilities. Namely, the process stays in the state j for some random time interval Tj (see 
Fig. [TJ distributed in accordance with the probability density 4>j(r) and jumps at the end 
of each residence time interval into another state i with the probability Py (Pjj = 0). The 
sum of these probabilities must add to one, X^i-^j = 1- Moreover, the time- intervals are 
assumed be independently distributed. Then, the just given trajectory description (Fig. [TJ 
defines the considered non-Markovian process completely. Furthermore, one assumes that 
for each given conformation the corresponding ET rate is constant, kj := 

The task is to average the relaxation dynamics in Eq. ((TJ over the stationary confor- 
mational fluctuations of the ET complex assuming an initial preparation of the electronic 
subsystem always in the excited donor state, i.e., Pd(0) = 1. For the stochastic process 
£(t) prepared at to — in a certain nonequilibrium state j a similar problem was solved in 
Refs . 2 ^ 21 . However, to arrive at the stationary averaging one must perform also an averag- 
ing over the stationary initial preparations of For non-Markovian processes this task 
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is not trivia l 17 i 22 i 23 i 24 . Indeed, by preparing the electronic subsystem at to = in one and 
the same electronic state with the probability one, Pd(0) = 1, the whole ET complex will 
be found in a state j with the stationary probability Pj(0) = pf = hm.t-^ 00 Pj(t). However, 
this state will be already occupied before t for some unknown and random time r*. The 
conditioning on and averaging over this time must be done in a stationary setup. This fact 
is well known in the theory of renewal processes 17 . The averaging over r* yields a different 
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RTD for the first time sojourn, ipj (t) = ipj(t)dt/ (tj) 17 , where (tj) = J °° T-ipj(r)dr 
the mean residence time (MRT) in the £j conformation. These MRTs are assumed to be 
finite in the following. Only for the strictly exponential RTDs (Markovian case), the RTD 
for the first sojourn ipj(r) is the same as ipj(r), i.e., ipf\r) = ipj(r). 

A similar problem of stationary averaging was solved for a stochastic Kubo oscillator^, 
x(t) = iou[£(t)]x(t), in Ref.—. By observing that the relaxation problem is mathematically 
equivalent to that of the Kubo oscillator with the imaginary frequency u — > ik, the use of 
the results from Ref.— yields for the (Laplace-transformed) stochastically averaged proba- 
bility (pD(s)) st = J °° (pD(t)) st exp(—st)dt to occupy the excited donor state the following 
expression: 

Pf v-^ 1 - ipj(s + kj) pf 



j 



l_^ (s + fc; ) / i \ 
1 -ip n (s + k n ) p% 



x P, 



S + k n \1~n) 



where the stationary occupational probabilities pf of the conformational states £j are the 
solutions of 

In Eq. ©, D nm (s) = 5 nm ip m {s + k m ), I is the unity matrix and f(s) := e~ st f(t)dt 
denotes the Laplace-transform 26 of any corresponding function f(t). 

The result in Eq. (J2J) provides the formally exact solution of the studied problem. The 
first sum in Eq. (j2J presents the solution of the problem with frozen static disorder. It can 
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be trivially inverted to the time-domain: 

TV 

(Pz?0)) stauc = J^pf exp(-%£) . (4) 

The remaining part presents nontrivial dynamical contributions. Its analytical inversion 
to the time domain is possible in a few exceptional cases only. However, in many cases it 
can reliable be inverted numerically by use of the Stehfest method^!. Moreover, the exact 
analytical expression for the mean survival (transfer) time^, tet '■= Jq° {PD{t))dt, follows 
as tet = (pD(fy)st- It can be readily given in an explicit form for ergodic processes with 
uniform mixing, where = (tj) / Xlj( r i) anc ^ T^j^ij = 1~- This is also the case of two- 
state conformational fluctuations with P\ 2 = P 2 \ = 1 (alternating renewal process^) which 
is investigated below. 



III. TWO-STATE CONFORMATIONAL FLUCTUATIONS 



In this simplest case which is of major interest in applications, the general result in Eq. 
(J2J) reduces to: 

(Pn(s))st= , * v jE-TF (5) 
(n) + (r 2 > I tf 9 s + k : 



3=1,2 



'j 



(hi - k 2 f [1 -Ms + h)][l - Ms + k 2 )\ 



}■ 



(s + hyis + k 2 y 1 _ + k x )Ms + h) 

This result can also be obtained from one in Ref.— by virtue of the above formal analogy 
with the Kubo oscillator. The corresponding exact expression for the mean transfer time 
reads: 

1 f in) M (r a > 

Tet = , ; , ; ; 1 ~, h 



(ti) + (r 2 ) I h k 2 

1 iy [1- & (%)][! -yi 2 (fc 2 )] 1 

>2 fcj 1-^1(^)^2(^2) J ' 

This expression is valid for arbitrary RTDs with finite MRTs and presents one of important 
results. It defines also an effective ET rate as r : — Tet * 



A. Markovian fluctuations 



Let us first reproduce the known results for the case of Markovian fluctuations from the 
above general formal expressions. For a Markovian process, the both RTDs are strictly 
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exponential: ifrxpfa) = (ti^) -1 exp(— t/(ti j2 )). The trajectory description with such RTDs 
is equivalent to a rate description with the conformational transition rates 71,2 = (71,2) ~ 1 . 
The corresponding Laplace-transformed RTDs read ipi,2 (s) = 1/(1 + s (t- 12 ) ) • In this simplest 
case, the inversion of Eq. (j5J) to the time domain yields a biexponential relaxation, 

(Pd (*))«* = ci exp(-rit) + c 2 exp(-r 2 t) , (7) 
with the relaxation (transfer) rates 

IY2 = ^{ki + k 2 + 71 + 72 =F VWi + 7i - h - 72) 2 + 47i7 2 } (8) 
and the weight coefficients (probabilities of different relaxation channels) 

Cl2 = j / 1=| _ 7i+72 + (fci-fc2)(7i-72)/(7i + 72) \ ^ 
2 I a/(A;i + 7i - h - 72) 2 + 47172 J 

The corresponding mean transfer time reads 

r = ^171 + fc 2 7 2 + (71 + 72) 2 , . 

ET (71 + 72) (hfa + &172 + fc 2 7i) ' 
The result in Eqs. (0)-© reproduces, e.g., the known result in Ref.— for the case of identical 
transition rates 71 = 72 = 7- This latter one can also be obtained from the well-known 
solution of the Kubo oscillato r 2 ^ 31 - 1 by virtue of the above mentioned formal analogy. The 
remarkably simple result in Eqs. (JZj)-© captures in essence three different transport regimes 
with dynamical disorder. If the conformational dynamics is slow, 7172 <C {k\ + 71 — k 2 — 
72) 2 /4, then the limit in Eq. (jlj) is approximately reproduced with pf 2 = (ti j2 )/ (fa) + fa)) 
(the regime of quasi-static disorder). If 712 ^> k\^ (the fast modulation regime), then ET 
is well described by the averaged rate (k) = pfki +p|*/c 2 , i-e., 

(p D (t)) st ^exp(-(k)t). (11) 

One more interesting regime occurs for 71 = 72 = 7 when jr^r «C 7 <C (k)/2. In this case, 
the mean transfer time tet equals approximately to the autocorrelation time r corr = 1/(27) 
of the Markovian two-state fluctuations, as it can be deduced from Eq. (|10J) . i.e. 

tet ~ T corr . (12) 

In other words, ET becomes gated, i.e., it is chiefly controlled by conformational 
fluctuation o 8 ! 10 ! 31 . An interesting gated regime emerges in the case where both ET rates 



are very different, k 2 <C ki, and the ET complex stays most of the time in the weakly 
conducting conformation, i.e., (r 2 ) ^> (n). Then, for k 2 <C 72 <C 71 <C /ci Eq. (fTUj) yields 



T E T ~ (T2) . 



(13) 



At the first look, the result in Eq. (|T3*j) might seem puzzling: the mean transfer time is given 
by the mean time intervals spent in the weakly conducting conformation. It is, however, 
quite natural since 72 = (72) _1 is nothing but the switching-on rate which corresponds to 
the rate limiting step of the overall ET process. 

B. Non-Markovian fluctuations 

To demonstrate a new surprising effect due to a broad non-exponential distribution of the 
residence times, let us consider the following tractable model. The conformation with the 
larger ET rate k\ (the "on" -conformation) has an exponential distribution of residence times 
with MRT (ri). In contrast, the "off" -conformation with a much smaller rate k 2 <C k\ is sup- 
posed to consist of a large number of quasi-degenerate sub conformations. Correspondingly, 
it can undergo a (sub) diffusion among these subconformations, which all possess one and the 
same ET rate k 2 , with a characteristic diffusional time r^. This time can be related to the 
volume of the corresponding conformational subspace (modeled here as an one-dimensional 
box, cf. Appendix EI), or to the number of quasi-degenerate diffusional substates, i.e. to 
the statistical weight of the conformation "2" . For this reason, within a statistical thermo- 
dynamic interpretation this time can also be related to the entropy of the conformational 
macrostate "2". Furthermore, one assumes that a conformational rate transition can occur 
into the "on" -conformation by crossing one of the boundaries of the underlying quasi-one- 
dimensional conformational subspace. An effective rate of such transition is 72 := (t 2 )~ 1 - 
The corresponding diffusional model, which is outlined in more detail in the Appendix El is 
characterized by the Laplace-transformed RTD 32 ' 34 
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(14) 



where 




(15) 
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The function g a (z) one has the following remarkable properties: g a (0) = 1 and g a {z) de- 
cays monotonously to zero with increasing z. The case a = 1 corresponds to the normal 
diffusion- 2 and < a < 1 corresponds to a subdiffusion with index The full details of 
the derivation of Eq. (J 14)) and the relevant discussion are given in Ref.— . 

Let us consider first the case of normal diffusion, a = 1. For r d = 0, Eq. (jT4"jl yields 
a single-exponential RTD. For T d <C (T2), the intra-conformational diffusion effects are not 
essential. However, for T d (r 2 ) the situation changes. Namely, for (r 2 ) 2 /T d <C r <C 
the corresponding RTD becomes a negative power law, ip 2 (r) oc t~ 3 / 2 , which ends by an 
exponential tail for r > r a 32 . The used model allows one to trace the influence of non- 
Markovian effects on the averaged transport process by assuming different variable values of 
T d with other model parameters, ki, k 2 , (n), (r 2 ), being fixed and defining the corresponding 
Markovian reference model. For this goal, we consider below a transfer regime with fci(ri) <C 
I5 k 2 (T 2 ) <C 1 which must be regarded as a fast modulation regime in the Markovian case. 
For fei(ri) <C 1, Eq. (JBJ) yields the following approximate analytical expression for tet- 

_ 1 A / \ , / / \ - k 2) 2 (ri}(T 2 }g a (k 2 T d )\ 

kik 2 {{Ti) + (r 2 » V h^) + k 2 {T 2 )g a (k 2 T d ) J 

Interestingly enough, for (n) = (t 2 ) = (r), the dependence on (r) in Eq. (fTBj) drops out 
and it becomes 

_ 1 f, , , (fci -fc 2 )V(fc 2 r d ) ^ 

r^T ~ 77] r~ Ui + fc 2 : -: f, r- • (17) 

2k x k 2 \ fci + k 2 g a {k 2 T d ) J 

The result in Eq. ()16p allows one to trace analytically the emergence of the regime of quasi- 
static disorder with the increase of r d . Namely, for r d = 0, g a (0) = 1 and Eq. (fTBj) yields 
tet ~ (k) -1 which corresponds to the averaged rate description in accordance with the 
intuitive expectations. However, if T d — ► 00, then g a {k 2 T d ) — > and r ST = pf/cf 1 +^2*^2- 
This corresponds to the quasi-static regime. Therefore, a transition to the quasi-static 
regime can occur due to an increase of the conformational diffusion time T d even if the 
MRTs (ti) and (r 2 ) are kept fixed well within the Markovian fast modulation regime. The 
corresponding criterion reads: 

g a ik 2 T d ) < 1. (18) 

This surprising result is corroborated by the numerical inversion of Eq. (jSJ) depicted in Fig. 
12] for the following set of parameters: k\ = 1, k 2 = 0.1, (ti) = (t 2 ) = 0.1 with varying r d . 
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The numerical inversion of the Laplace-transform (JHJ) was done using the improved Stehfest 



in the corresponding asymptotic series expansion 27 (N = 16 in the original and standard 
version) should correspond to the digital precision of calculations 28 . The latter one (along 
with the number of terms in the Stehfest series) was made variable (increased) until the con- 
vergency of the numerical results was reached. This was made possible by implementing the 
Stehfest algorithm using the Computer Algebra System MAPLE which allows for numerical 
calculations with an arbitrary precision. N = 32 was sufficient for the numerical calculations 
done in this paper to arrive at the required accuracy (within the line width in the figures 
representing the numerical results). Note that the mean residence times (t 1j2 ), which are 
used for calculations in Fig. |21 are indeed much less than the inverse of the corresponding 
ET rates. This relates also to the mean turnover time, Tt U m '■= fa) + fa) <C fcj" ,A;J . 
Therefore, one intuitively expects that many stochastic turnovers should occur before the 
electron donor population changes significantly. This in turn should validate an average 
rate description of the ET process. Contrary to this intuitively clear picture which is well 
established in the case of Markovian fluctuations^, an increase of results in a gradual 
transition to the regime of quasi-static disorder, see in Fig. |21 This presents certainly 
a paradox: Judging from the time scale r iurn the ET complex fluctuates rapidly (on the 
time scale of decay of the donor population) between two conformations. Nevertheless, the 
averaged rate description breaks down and for large the averaged ET process becomes 
non-exponential entering a slow modulation regime. It can be even better described by the 
quasi-static approximation (top curve in Fig. EJ) which corresponds to a "frozen" disorder. 

A deceivingly clear resolution of this paradox is based on the observation that for the con- 
sidered non-Markovian process the effective autocorrelation time can be much larger than 
the autocorrelation time of its Markovian counterpart with the same MRTs^ 3 .. The mean 
autocorrelation time r COTr can serve as a quantifier for memory effects. It can be defined 
as r corr : = J" °° \cov(t)\dt, where cov(t) := (5k(t)5k(0)) st /(8k 2 ) st is the normalized autocor- 
relation function of k(t) (here, 5k(t) := k(t) — (k)) and obtained by using the Stratonovich 



for the Laplace-transform of the normalized autocorrelation function of a stationary two- 



metho d 27 ' 28 following a "rule of thumb" found empirically: the number N of terms used 



formul a 14 ' 33 ' 35 




(19) 
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state renewal process. If cov(t) does not change its sign in time (the present situation), then 
Tcorr = lim s ^ cov (s) and by use of Eq. [191 the mean autocorrelation time reads— 

Tcorr = ^(^l + ^i) T M (20) 



where 



G> = (21) 

is the coefficient of variation of the corresponding RTD ipj(r) and 

(ti) + (r 2 ) 

is the autocorrelation time of the Markovian process with the same MRTs (ti) and fa). 
For the considered model, r corr = (1 + Td/Sfa))TM- It can be much larger than tm with the 
increase of r d . In contrast with the Markovian process, the time-correlations do not decay 
on the time scale of stochastic turnovers and can persist very long. This can be considered 
as a reason for the transition to the quasi-static regime. Judging from the autocorrelation 
time scale it can indeed be the quasi-static regime. 

Accepting this reasoning uncritically, one might conclude that for the processes with a 
divergent RTD variance, C2 —>■ 00, and thus with r corr —>■ 00, the quasi-static regime will be 
realized always. However, such a conclusion is not correct. Indeed, let us take some a < 1 
in Eq. (fTlj) (this corresponds to a subdiffusion with index oM-). Then, it is easy to show 
that formally r corr = 00 for all 7^ and one might naively expect that the regime of 
quasi-static disorder in the ET kinetics will readily be realized. However, the numerics do 
not support this expectation. For example, for a = 0.1 and for = 100 the averaged ET 
is not exponential indeed. Nevertheless, it is closer to the regime of fast modulation rather 
than to the regime of quasi-static disorder. This non-occurrence of the regime of quasi-static 
disorder agrees, however, with the criterion (|18|) which yields in this case g a (k2Td) ~ 0.72 
for the given parameters. This value is indeed closer to g a {k2Td) = 1 (the averaged rate 
description) rather than to g a {k 2 T d ) = (the regime of quasi-static disorder). 

It is not the divergence of the mean autocorrelation time which does play a crucial role, 
but an intermittent, bursting character of non-Markovian fluctuations which is important for 
the discussed phenomenon to occur. For r d = 100 in Fig. the conformational fluctuations 
occur by strongly correlated bursts separated by long interburst periods as revealed by the 
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stochastic simulations of k(t) (see in Fig. Efa),(b) where this striking feature is presented 
for a different model). Namely, even if the mean residence time (r 2 ) in the weakly conduct- 
ing conformation is very short, (r 2 ) <C k^ 1 , the ET complex can nevertheless be trapped 
occasionally in this conformation for a very long time period of the order of ^> k^ 1 . The 
time Td is a characteristic time which is required to explore by random walk the manifold 
of quasi-degenerate conformational substates within the corresponding macroconformation. 
Such long time intervals separate the bursts within which the ET complex fluctuates rapidly 
between its two macroconformations. The mean residence time (r^) of the time intervals 
spent in the conformation £2 (conformation "2") within a burst could obviously be much 
less than the overall (r 2 ) (since this latter one accounts also for the long interburst sojourn 
in the same conformation), i.e. (r^) <C (t 2 ). For (ti) = (r 2 ) and pf = p 2 4 , this means 
that during the the burst duration the ET complex spends most of the time in the confor- 
mation "1" since (ti) 3> {t' 2 ). If the burst durations (r burst ) are, on average, much larger 
than the inverse ET rate k± , i.e. (r^rst) 3> k± , then it becomes clear why the quasi- 
static averaging does apply in such a situation, even if (ti) = (r 2 ) <C , ftj 1 . This is the 
reason why the averaging over many repeated realizations of electron transfer in the same 
ET complex, or an ensemble averaging for a number of such complexes done at the same 
time can indeed yield a quasi-static regime even if a large portion of these complexes does 
fluctuate very fast on the time scale of ET while another one stays temporally "frozen" in 
the "off" -conformation. However, even if r corr is formally infinite (a < 1), the fluctuations 
do not necessarily have a distinct bursting character. Then, the discussed phenomenon does 
not necessarily occur. This qualitative explanation is fully confirmed within the following 
exactly solvable non-Markovian model. 

C. Exactly solvable non-Markovian model 

The simplest non-Markovian model is obtained by assuming a biexponential RTD, 

^ 2 (r) = Octi exp(-«ir) + (1 - 6)a 2 exp(-a 2 r), (23) 

where «i j2 are two different transition rates for the conformational transition from the con- 
formation "2" to the conformation "1" . They can be realized with the probabilities < 9 < 1 
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and 1 — 9, correspondingly. The corresponding mean residence time is 



(r 2 ) = 9/a 1 + (l-9)/a 2 , 



(24) 



and the coefficient of variation ()21|) for this distribution is 




(25) 



where the functions fi t2 (9,() are given in Eq. (jB6|) and ( = a\ja 2 is the ratio of two rate 
constants in Eq. (|2*H|) . For ( <C 1 and ( ^ 9 the following approximate expression is valid 



Note that C 2 can be very large. This applies when 9 <C 1 and £ <C 1, see in Fig. El This 
is the case of distributions (J2*3|) possessing a very small, but very broad tail. Furthermore, 
the Laplace-transform in Eq. (J19)) can be inverted analytically for this model to obtain a 
biexponentially decaying autocorrelation function cov(t). In the parameter regime corre- 
sponding to Eq. ()26|) . this autocorrelation function becomes, however, approximately single 
exponential, 



with tm in Eq. (}2*2*j) . The result in Eq. ()27|) becomes exact in the formal limit a 2 — > oo, 

where if) 2 {r) — > 9ct\ exp(— air) + (1 — 9)5(r) and (t 2 ) — > 9j<X\. This striking feature must be 
emphasized: the exists a whole family of non-Markovian two-state processes with an expo- 
nentially decaying autocorrelation function. The autocorrelation time of these processes is, 
however, enlarged by the factor 1/9 as compare with their Markovian counterpart possessing 
the same MRTs (ti) and (r 2 ). In accordance with Eq. (|2*7|) . the autocorrelation time of such 
non-Markovian processes tm/9 can be much larger than the autocorrelation time tm of their 
Markovian counterpart. 

The considered two-state non-Markovian process can also be embedded into a three state 
Markovian kinetic scheme with two substates (subconformations) corresponding to the con- 
formation "2". The two-state fluctuations of the ET rate k(t) present nevertheless a non- 
Markovian process. The just described model can be solved exactly. Namely, Eq. (jSJ) yields 
after some lengthy calculations a rational function 




(26) 



cov(t) ps exp(-8t/T M ) 



(27) 



(Pd(s)) 



s 2 + a%s + a 



(28) 



s 3 + b 2 s 2 + bis + &o' 
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where the polynomial coefficients dj, hi are given in the Appendix[Bl The analytical inversion 
of Eq. ()28|) to the time domain yields formally a three exponential kinetics 

(p D (t)) st = 

Ciexp(-Tit) (29) 

i=l,2,3 

with the rates Tj which are given in the Appendix |Bl Eqs. ()B7|) - (jB9|) . and the weights 
Cj = (rf — aiTi + ao)/(3rf — 2b 2 Ti + b\). A numerical evaluation for all cases studied in 
this paper shows, however, that the ET kinetics remains effectively biexponential (at most) 
since one of the rates T 2 exceeds much the other two and the corresponding weight c 2 is 
negligible small. 

In order to establish a criterion similar to Eq. (fTHjl . the Laplace-transform of (}2"3*|) can be 
represented in a form similar to Eq. ()14|) . 

1 + s(T2)g{s{T2),9,Q 

where 

q(z ec)= [0 + (i-e)c? + Cz m 

yK ' ' w [# + (i-#)C] 2 + [# + (i-#)C][i-# + flCK 

and MRT (r 2 ) is given in Eq. (J24|) . Then, the criterion (|18|) is get modified for the present 
model as 

<7(M*2>,0,C)<1 (32) 

with ^2(72) *C 1. Let us clarify if it is possible to fulfill the criterion (|32|). It is rather 
nontrivial in the present case since we do not have such an independent parameter as the 
diffusional time which can easily be adjusted to satisfy the condition (|18|) for a sufficiently 
large value of r^. Therefore, a deliberate choice of parameters is required. We consider the 
following set of parameters k± = 1, k 2 = 0.2, and (71) = (r 2 ) = 0.2 to be still within a 
fast modulation regime in accordance with the Markovian criteria. The plot of the function 
g (z, 6, C) versus 9 at z — 0.04 and at the different values of ( in Fig. |1] reveals that the 
criterion (J32)) can indeed be fullfield for a very small ratio of the rate constants, starting 
from £ = 10~ 3 , and for a very small, but finite splitting probability 9 (note the presence 
of an optimal value of 9 in Fig. IH where g (fc 2 (r 2 ), 9, () reaches a conditional minimum). 
The corresponding conformational dynamics is profoundly bursting indeed, with long time 
intervals between bursts as it is revealed by the stochastic simulations shown in the Fig. El 
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(a),(b). Here, the interburst intervals have the characteristic time scale Ti n t er ~ ( r 2)/Q = 
20 > ki2- The burst durations have about the same characteristic time scale and within 
the burst the excursions into the state "2" are extremally brief and the process stays mostly 
in the state "1". These are the reasons why a regime similar to the regime of quasi-static 
disorder is realized. In the case 9 = 0.5 (which is close to the regime of fast modulation, cf. 
Fig. H)J), k(t) has no such profoundly bursting character - compare the Fig. 5(c), where a 
sample trajectory of this non-Markovian process is depicted, with the Fig. 5(a), (b) and with 
the Fig. 5(d), where a sample trajectory of the corresponding Markovian process is shown. 
Clearly, in this case k(t) reminds strongly its Markovian counterpart. This similarity can 
be, however, somewhat deceiving. The sample trajectories in Fig. 5(c) and Fig. 5(d) are 
replotted in Fig. 5(e) and Fig. 5(f), respectively, on a much shorter time scale. The both 
processes have the same (n) = (r 2 ) = 0.2. However, the trajectory in Fig. 5(e) corresponds 
to a non-Markovian process, where some very brief excursions occurs into the state "2" 
during the sojourn in the state "1". This non-Markovian process has, in accordance with 
Eq. (J2*7|) . the autocorrelation time which is twice bigger as compare with its Markovian 
counterpart in Fig. 5(f). 

The stationary stochastic trajectories in Fig. El were generated due to the following 
algorithm. First, a random number < y\ < 1 is drawn from the uniform distribution^ 
and one decides in which conformation to start. Namely, if y± < pf, then one starts in the 
state "1", otherwise the state "2" is taken as initial. If the initial value is k\, then a random 
residence time-interval T\ is drawn from the corresponding exponential distribution. Here, 
in accordance with the rule of transformation of random variables, T\ = (n) \n(l/y), where 
y is a new random number taken from the uniform distribution 3 -^. After this, the process 
k(t) flips into the state "2" and two new random numbers y\ and y2 are drawn from the 
uniform distribution. y± is used to decide in accordance with which exponential distribution 
■021 ( r ) = Q!i exp(— Qir), or ip-nij) = a2exp(— a 2 r) the process will sojourn in the state "2". 
If Hi < 9, then it resides in the state "2" in accordance with the distribution ip2i{T~), i-e., 
r 2 = ln(l/y 2 )/«i, otherwise the distribution ?/> 22 (t) is used, i.e., r 2 = ln(l/?/ 2 )/a; 2 . After 
the sojourn of duration r 2 in the state "2", the process flips into the state "1" and a new 
random number y is generated to obtain the residence interval T\ = (ti) ln(l/y). The whole 
procedure repeats. If the initial state happens to be the state "2" (cf. Fig. 5(b)), then the 
whole procedure is modified. Namely, for the first sojourn in the state "2", a different RTD 
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( r ) = /°° ip2(t)dt/ (r 2 ) is used (see the discussion in Sec. |TTJ) which is 



^2 CO = @i a i exp(— air) + (1 - #i)a 2 exp(-a 2 r) 



(33) 



with 











(34) 



+ (i-0)C 



This modification is very important. The generated stochastic trajectories will not be sta- 
tionary otherwise. As a matter of fact, Q\ can be even closer to one when ^ < 1, but 
C <C 9. The stationary probability to find the considered ET complex initially in the "off "- 
sub conformation with the transition rate ol\ <C a 2 is and not p?f9\ 9 is in fact the 

splitting probability, i.e. the probability to switch from the macrostate "1" into the sub- 
state of the macrostate "2" with the conformation transition rate ai. The disregarding of 
this very important point would result in a serious mistake as the generated trajectories 
would not be stationar y 17 ' 22 ' 24 and thus would not correspond to the stationary conforma- 
tional fluctuations which are assumed throughout this work. The issue of noise stationarity 
is crucial for our consideration. Similar to the well-known Gillespie algorithm 4 ^, the used 
here stochastic algorithm to generate the trajectory realizations of the considered stochastic 
process is exact. 

Generally, the minimum of g (fc 2 (r 2 ), 9, () is not related to the maximum of C 2 (0, C) (and 
thus to the maximum of the mean autocorrelation time r corr ), even if both extrema take 
place in the same parameter region <C 1,( < 1, compare Fig. Eland Fig. 0J Therefore, 
the proper criterion is not related to the mean autocorrelation time of the considered non- 
Markovian process - this being in accordance with the discussion in the subsection IHI Bl 
The validity of the criterion ()32j) is fully confirmed in Fig. El for ( = 10 -3 and three different 
values of 9. Namely, for 9 = 0.5 the ET kinetics is close the regime of fast modulation, 
whereas for 9 = 0.01 it is close to the regime of static disorder. With a further decrease 
of 9 the kinetics starts actually to depart from its quasi-static limit (see the dotted line 
for 9 = 0.001 in Fig. El even if the coefficient of variation C 2 and, therefore, the mean 
autocorrelation time Td attain for this value of 9 their maxima, cf. Fig. El). 
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IV. SUMMARY AND CONCLUSIONS 



Let us summarize the main results of this study. A simple model of rate processes with 
non-Markovian dynamical disorder has been investigated where non-Markovian rate fluc- 
tuations were modeled by a renewal process of the continuous time random walk type. A 
special attention has been paid to finding the averaged dynamics with averaging done over 
the stationary noise realizations. Such averaging corresponds to a stationary environment. 
For the considered class of non-Markovian processes the main result is given in Eq. (j2J). 
In the simplest case of two-state rate fluctuations it reduces to the results in Eqs. (0) and 
(JBJ). All these novel results are exact. The previously known results for the case of two-state 
Markovian fluctuations are reproduced in the simplest particular case which corresponds 
to the exponentially distributed residence times in the trajectory description of the consid- 
ered stochastic process. This provides a validity test for the developed theory. The general 
theory was applied to study two important non-Markovian models: the model incorporat- 
ing a conformational diffusion within one of the macro conformations over the energetically 
quasi-degenerate microstates and the model assuming that one of the macroconformations 
consists of two subconformations only (the simplest, minimal non-Markovian model). The 
former model yields a power law distribution of the residence times in the corresponding 
macroconformation and the latter one yields a biexponential distribution of the residence 
times. 

Within the model featuring a power law distribution, a striking paradoxical effect was 
revealed. Namely, the transport process can approximately be described by the quasi-static 
approximation notwithstanding that its rate fluctuates rapidly on the transport time scale 
judging from the mean residence time intervals spent by the charge transferring macro- 
molecule in its corresponding macroconformations. A proper explanation of this paradoxical 
effect is given in terms of the bursting features of the underlying non-Markovian fluctuations. 
This explanation was confirmed and reinforced within the latter, simplest non-Markovian 
model which was solved exactly and the corresponding parameter regime was found where 
the discussed paradoxical effect exists as well. The proper analytical criterion for its oc- 
currence was obtained for both models in similar analytical forms, cf. Eq. ()18j) and Eq. 

In the same spirit, the corresponding criterion can be found for any other model of the 
residence time distribution ^(t) with a finite mean residence time fa). This is so because 
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any such distribution (more precisely, its characteristic function, or Laplace-transform) can 
be written in the form of Eq. (JT3j) with the corresponding function g(s, {A}) found, which 
depends on the set {A} of the model parameters such that g(0, {A}) = 1. If there exists such 
a particular parameter set (A)} that g(k 2 , {^4o}) <C 1 for /c2( r 2) <C 1 and A; 1 (r 1 ) <C 1, then 
our theory predicts the occurrence of the discussed paradox of non-Markovian dynamical 
disorder for this particular parameter set, cf. Eq. (JTfij) and the corresponding discussion 
below Eq. (fTHjl . 

In practice, the paradoxical effect revealed in this work implies that if an ensemble mea- 
surement indicates that a charge transfer (or another relaxation process) takes place in an 
apparently quasi-static, biexponential regime this does not mean yet that we are dealing 
indeed with an ensemble of ET complexes "frozen" randomly in one of the two different 
conformations. The following situation is possible: one-half of all ET complexes is indeed 
"frozen" temporally in one of the conformations while another half (the proportions here can 
be different) fluctuates very rapidly between different conformations. In such a situation, 
the repeated measurements only done with single molecules can help to uncover their actual 
conformational dynamics. This would require but a sufficient time resolution of the mea- 
surement apparatus since otherwise the most of switching effects can be missed merging in 
an apparently long sojourn without resolving its actual fine structure. 

In conclusion, our study reveals that the non-Markovian effects can have a profound im- 
pact on the rate processes with dynamical disorder. Their role is expected to be essential for 
the charge transfer processes in macromolecular complexes with an intrinsic conformational 
dynamics. 
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APPENDIX A: A MODEL OF CONFORMATIONAL (SUB)DIFFUSION 



The conformational subspace x of the state "2" is modeled by an one-dimensional box of 
the length L, —L < x < 0. The conformational diffusion within this box is assumed to be 
described the fractional diffusion equation (see, e.g., a review 3 ^ and a popular account^ 7 -, as 
well as the further references therein). The fractional diffusion equation can be written in 
the form of standard continuity equation for the probability density P(x,t), 

dP(x,t) dJ(x,t) 

with the probability flux J(x, t) being modified due to the fractional time-derivative, 

w x t j*~ a1-« dP(x,t') , . nS 

J{x, t) = -K a D t ^ ; . (A2) 

In flA2|) . Q D t (...) = fj^jM Jo (t-vy-<* * s ^ ne integro-differential operator of fractional 
derivative introduced by Riemann and Liouvill o 36 ^ 37 and K a is the fractional diffusion co- 
efficient. For free space, — oo < x < oo, the above fractional diffusion equation describes 
an anomalous subdiffusion, (x 2 (t)) oc t a , with the index of subdiffusion < a < l 36 ' 37 . 
For a = 1, it reduces to the conventional diffusion equation. In the present case of the 
conformation diffusion restricted in space, the diffusional time Td '■— (L 2 / K a ) l l a becomes an 
important model parameter. 

Furthermore, the left boundary x = —L is assumed to be reflecting, J(—L,t) = 0, 
while the right boundary is radiative, J(0,t) = LP(0,t)/ (r 2 ) (here the transition into the 
state "1" occurs with the effective rate 72 = (t2) _1 ). Since the backward transition from 
the state "1" to the state "2" occurs by crossing this right boundary in the back direc- 
tion, to obtain the corresponding residence time distribution ^(t) one must solve first 
the fractional diffusion equation with the stated boundary conditions complemented by the 
initial condition P(x, 0) = 5(x — Xq) with xq — > 0_. Then, the RTD ip2{ T ) follows as 
4>2( T ) — —-j^ f^. L P( x , T )dx. This task has been accomplished in Ref. 34 using the Laplace- 
transform method. The result is given in Eq. (fT^j) . An analogous problem in the case of 
normal diffusion was solved in Ref.- 3 ^. 

The used here continuous diffusion model has been derived in Ref.^ from a CTRW gen- 
eralization of the discrete state diffusion model originally proposed by Millhauser, Salpeter, 
and Oswald 33 to explain some features of the conformational gating dynamics of biological 

18 



ion channels. It can serve, however, to model the conformational diffusion in other proteins 
as well. 



APPENDIX B: POLYNOMIAL COEFFICIENTS AND THE RATES I\ 

In this Appendix, the exact solution of the model of subsection IHI Cl is presented. Namely, 
the polynomial coefficients in Eq. (J28|) are: 

, , h(r 2 } + h(n} l l l . f . 

ai = k 2 + , , , . + t — r + 7 — r + 7— r/i(0, C , BI 
(n) + (t 2 ) (ri) (t 2 ) (t 2 ) 

. ^1(7-2) + fc 2 (n) 1 1 

a = k 2 — / — v / — v V 1 — r + 



(n) + (r 2 ) (n> (r 2 ) 
+ A ^ + "^ WX) + A f A + A ) JWC), (B2) 



(7-2) (r 1 ) + (r 2 ) (r 2 ) V(n) (r 2 ) 

6 2 = 2^ + ^ + A + A + A^'O. (B3) 
(ti) (t 2 ) (t 2 ) 

&i = fc 2 + 2A; 2 (fc + A) + A^ + k ^ ^ + /i(*'0) ( B4 ) 
V in) J (r 2 > 

1 

+ -r- 



>+ 




)/a' 


1 

(ri> 


) + 


1 

(7-2) 


) + 


(t 2 ). 





6 = A; 2 ( fci 

+ A (-TT + 7^)^.0- (Br,) 



In Eqs. (jBl|l -flB jl . two auxiliary functions, 

h{ex) = (i-ex + e -, 

W,0 = ({i-e)y/c + ^j , (B6) 

are introduced. Furthermore, ( := a\ja 2 is the ratio of two "on"-rates, corresponding to 
the two sub conformations within the macro conformation "2" and the mean residence time 
(t 2 ) is given in Eq. 

The rates Tj in Eq. (|29jl are given by the roots (with the change of the sign) of the cubic 
equation which is obtained by setting the denominator in Eq. ()28|) equal to zero. These 
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rates read explicitly, by use of the well-known trigonometric solution of the cubic equation^ 6 .: 

(B7) 









/ Q cos 


r\ = 


\ 

3 


- 2 \ 


r 2 = 


\ 

3 


- 2\ 


/ Q cos 


r 3 = 


\ 

3 


- 2 \ 


/ Q cos 



V 
3 

77 + 2n 
3 

1] + 4-7T 



(B8) 
(B9) 



where 



77 = arccos > ( B1 °) 

Q = (6 2 /3) 2 - &i/3, (Bll) 

R = (6 2 /3) 3 - hh/6 + 6 /2, (B12) 

and one assumes Q > and < Q 3 ^ 2 - 
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FIG. 1: A sample trajectory of the considered non-Markovian process. When the ET complex 
is prepared in the excited donor state D* at to = 0, the initial conformation j (here j = 1) was 
already occupied before to f° r some unknown time tJ. The conditioning on and the averaging 
over r* yield generally a different RTD ipj (r) of the first sojourn in the corresponding initial 
conformation. For the strictly exponential ipj{ T ) only (Markovian case), ipj '(r) = ipj( T )- 

FIG. 2: The stochastically averaged decay of the donor population {pD(t))st is plotted versus 
the time t (in arbitrary units) for several different values of the conformational diffusion time t&. 
An increase of Td (from bottom to top) causes transition from the single-exponential ET regime 
of fast modulation (solid bold line, = 0) to a non-exponential ET regime of the quasi-static 
disorder (dashed bold line). The mean residence times in the conformations corresponding to the 
rate k\ = 1 and the rate &2 = 0.1 are (t\) = 0.1 and (T2) = 0.1, respectively. 

FIG. 3: The coefficient of variation C2 of the biexponential residence time distribution (|23|) is 
plotted versus 6 for different ratios £ = axjai of the rate constants a\ and 02- The dotted line 
corresponds to the approximation in Eq. (|26|) , 

FIG. 4: The criterion function g(z, 9, in Eq. ()31j) is plotted versus 9 for different ratios C = 01/02 
and for the fixed value z = 0.04. 

FIG. 5: Sample trajectories of the two-state ET rate fluctuations: (a)-(c) and (e) corresponds to 
non-Markovian rate fluctuations with the single exponential distribution of the residence times in 
the state "1" with MRT (n) = 0.2 and with the biexponential distribution ()23|) of the residence 
times in the state "2" with MRT (r 2 ) = 0.2; ( = 0.001 and 9 = 0.01 (in (a) and (b)) versus 9 = 0.5 
in (c) and (e). In (d) and (f) a sample realization of the Markovian process with the same MRTs 
is depicted on two different time scales. 
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FIG. 6: The stochastically averaged decay of the donor population {pD(t)) s t is plotted versus the 
time t (in arbitrary units) for the simplest non-Markovian model with the biexponential distribution 
(|23|). The solid bold line corresponds to the ET regime of fast modulation and the dashed bold line 
corresponds to the regime of quasti-static disorder. The actual ET kinetics lies in between these two 
limits. The following parameters are used in calculations: k\ = 1, A% = 0-2 and {t±) = fa) = 0.2, 
£ = 0.001. For 9 = 0.01 the ET kinetics is more close to the regime of quasi-static disorder and for 
6 = 0.5 the ET kinetics is more close to the fast modulation regime. 
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